CHAPTER 9 FIRST-ORDER DIFFERENTIAL EQUATIONS 



9.1 SOLUTIONS, SLOPE FIELDS AND EULER'S METHOD 

1. y' = x + y => slope of 0 for the line y = — x. 

For x, y>0, y' = x + y=> slope > 0 in Quadrant 1. 

For x, y<0, y' = x + y=> slope < 0 in Quadrant III. 

For |y| > |x|, y>0, x<0, y' = x + y=> slope > 0 in 
Quadrant II above y = — x. 

For |y| < |x|, y>0, x<0, y' = x + y=> slope < 0 in 
Quadrant II below y = — x. 

For |y| < |x|, x>0, y<0, y' = x + y => slope > 0 in 
Quadrant IV above y = — x. 

For |y| > |x|, x>0, y<0, y' = x + y=> slope < 0 in 
Quadrant IV below y = — x. 

All of the conditions are seen in slope field (d). 

2. y' = y + 1 => slope is constant for a given value of y, slope 
is 0 for y = — 1 , slope is positive for y > 1 and negative for 
y < — 1. These characteristics are evident in slope field (c). 



3. y' = — ^ =>• slope = 1 on y = — x and — 1 on y = x. 
y' = — * =>■ slope = 0 on the y-axis, excluding (0, 0), 
and is undefined on the x-axis. Slopes are positive for 
x > 0, y < 0 and x < 0, y > 0 (Quadrants II and IV), 
otherwise negative. Field (a) is consistent with these 
conditions. 



4. y' = y 2 — x 2 =>■ slope is 0 for y = x and for y = — x. 
For |y| > |x| slope is positive and for |y| < |x| slope is 
negative. Field (b) has these characteristics. 
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7. y = — 1 + / i X (t-y(t))dt^g=x-y(x);y(l) = -l + (t - y(t))dt = -1; g = x - y, y(l) = -1 

8 - y = /> ^ I = i; y(i) = ft = o ; g = 7 , y (t) = o 

9. y = 2- J o X (l + y(t))sintdt => g = -(1 + y(x))sinx; y(0) = 2 - J 0 °(l + y(t))sintdt = 2; g = -(1 +y)sinx, 
y (0) = 2 

10. y = 1 + / 0 X y(t) dt => g = y(x); y(0) = 1 + f° y(t) dt = 1; | = y, y(0) = 1 

11. Yi = Yo + (l - l) dx = -1 + (1 - =±) (.5) = -0.25, 

y 2 = yi + (l - dx = -0.25 + (l - ) (.5) = 0.3, 

y 3 = y 2 + (l - g) dx = 0.3 + (1 - f ) (.5) = 0.75; 

g + (i) y = 1 => P(x) = i , Q(x) = 1 => f P(x) dx = Jl dx = In |x| = In x, x > 0 =>• v(x) = e lnx = x 

=> y = 1 fx ■ 1 dx = i (f + c) ; x = 2, y =%1 =>• -1 = 1 + § => C = -4 => y = § - 7 

=S> y(3.5) = ^ - 33 = ^ ~ 0.6071 

12. yi = yo + x 0 (1 - y 0 ) dx = 0 + 1(1 - 0)(.2) = .2, 

y 2 = yi + xi (1 - yi) dx = .2 + 1.2(1 - ,2)(.2) = .392, 
y 3 = y 2 + x 2 (1 - y 2 ) dx = .392 + 1.4(1 - ,392)(.2) = .5622; 

= x dx => — ln|l — y| = y+C;x=l,y = 0 => -In 1 = 1 + C => C = - 7 => ln|l-y| = -^ + 7 
=> y = 1 - e (1 ~ x2)/2 =*► y(1.6) « .5416 

13- yi = y 0 + (2x 0 y 0 + 2y 0 ) dx = 3 + [2(0)(3) + 2(3)](.2) = 4.2, 

y 2 = Yl + (2x iyi + 2 yi ) dx = 4.2 + [2(.2)(4.2) + 2(4.2)](.2) = 6.216, 

y 3 = y 2 + (2x 2 y 2 + 2y 2 ) dx = 6.216 + [2(.4)(6.216) + 2(6.216)](.2) = 9.6969; 

g = 2y(x +1) => f = 2(x + 1) dx =>• In |y| = (x + l) 2 + C; x = 0, y = 3 ^ In 3 = 1 + C => C = In 3 - 1 
=> In y = (x + l) 2 + In 3 - 1 =*► y = e ( x+1 > a+ln3 - 1 = e ln3 e x2+2x = 3e x ( x+2) => y(.6) « 14.2765 

14. y x = y 0 + yg(l + 2x 0 ) dx = 1 + 1 2 [1 + 2(-l)](.5) = .5, 
y 2 = yi + y?(l + 2x0 dx = .5 + (,5) 2 [1 + 2(-.5)](.5) = .5, 
y 3 = y 2 + y|(l + 2x 2 ) dx = .5 + (,5) 2 [1 + 2(0)](.5) = .625; 

g = (1 +2x)dx =» - I = x + x 2 + C;x = -l,y = 1 => -1 = -1 +(-l) 2 + C => C = -l => 7 = 1 - x- 

^ y = i - x - x 2 ^ y(-^) = i — .5 — (. 5) 2 = 4 
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15. yi = y 0 + 2x 0 e x ° dx = 2 + 2(0)(.l) = 2, 

y 2 = Yl + 2xie x ‘ dx = 2 + 2C.l)e l2 (.l) = 2.0202, 
y 3 = y 2 + 2 x 2 e x 2 dx = 2.0202 + 2(.2) e 22 (.l) = 2.0618, 

dy = 2xe x2 dx => y = e x2 + C; y(0) = 2+>2=l+C+>C=l =>y = e x2 + l =+ y(.3) = c 3 ' + 1 « 2.0942 

16. yi = y 0 + (yo e x °) dx = 2 + (2 • e°) (.5) = 3, 

y 2 = yi + (yi e Xl ) dx = 3 + (3 • e 05 ) (.5) = 5.47308, 

y 3 = y 2 + (y 2 e X2 ) dx = 5.47308 + (5.47308 • e 10 ) (.5) = 12.9118, 

= y e x +> ^ = e x dx =>• ln|y| = e x + C; x = 0, y = 2 => In 2 = 1 + C => C = In 2 — 1 =>• ln|y| = e x + In 2 — 1 

=> y = 2e eX - 1 => y(l .5) =2e e ' 5 - 1 « 65.0292 



17. yi = 1 + 1 ( . 2 ) = 1.2, 

y 2 = 1.2 + (1.2)(.2)= 1.44, 
y 3 = 1.44 + (1.44)(.2) = 1.728, 
y 4 = 1.728 + (1.728)(.2) = 2.0736, 
y 5 = 2.0736 + (2.0736)(.2) = 2.48832; 

^ = dx => In y = x + Ci +> y = Ce x ; y(0) =1 => 1 = Ce° +> C = 1 => y = e x => y(l) = e « 2.7183 

18. yi = 2 + (j) (.2) = 2.4, 

Y 2 = 2.4 + (j=|) (.2) = 2.8, 
y 3 = 2.8 + (f|) (.2) = 3.2, 
y 4 = 3.2+ (f|)(.2) = 3.6, 
y 5 = 3.6 + (jj) (.2) = 4; 

f = f =>• In y = In x + C +> y = kx; y(l) = 2+>2 = k+>y = 2x+> y(2) = 4 



19. yi = -1 + 
y 2 = --5 + 



(-+ 

L vd J 

(— -5) 2 



vd^J 

y 3 = -.39794 + 
y 4 = -.34195 + 



05) = -.5, 

(.5) = -.39794, 

(.5) = -.34195, 
(.5) = -.30497, 



(— .39794) 2 
(— .34195) 2 



y/l3 

y 5 = -.27812, y 6 = -.25745, y 7 = -.24088, y 8 = -.2272; 

| = =► -±=2 > /i + C;y(l) = -l => 1=2 + C =+ C = 



1 ^ y_ 1-+A ^ y(5)_ l-2x/5 



.2880 



20. yi = 1 + (0 • sin 1 ) (i) = 1, 

y 2 = 1 + (1 -sinl) (i) = 1.09350, 
y 3 = 1.09350 + (| • sin 1.09350) Q) = 1.29089, 

y 4 = 1.29089 + (| • sin 1.29089) (§) = 1.61125, 

y 5 = 1.61125 + (| • sin 1.61125) Q) = 2.05533, 

y G = 2.05533 + (f • sin 2.05533) Q) = 2.54694; 

y' — xsiny +• cscydy = xdx => — ln|cscy + coty| = ^x 2 + C => cscy + coty = e - 5 x2+c = Ce^J x " 
=> 1 ^5° sy = Ce~i x " => cot( y ) = Ce _ J x2 ; y(0) = 1 => cot ( ) ) = Ce° = C =>■ cot ( 2 ) = cot(i)e _ i x2 

=> y = 2 cot -1 (cot(|)e _ * x '), y(2) = 2cot _I (cot(j)e~ 2 ) = 2.65591 
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21. y = -1 - X+ (1 + x 0 + yo)e x x ° => y(x 0 ) = -1 - x 0 + (1 + x 0 + y 0 )e x ° x ° = -1 - x 0 + (1 + x 0 + y 0 )(l) = yo 
| = -l + (l+xo + yo)e x - x »^y=-l-x+(l+x 0 + yo)e x - x <> = | - x =► g =x + y 

22. y' = f(x), y(x 0 ) = y 0 y = J x * f(t)dt + C, y(x 0 ) = J x *° f(t)dt + C = C=^C = y 0 =>y = £ f (t) dt + y 0 



23-34. Example CAS commands: 

Maple : 

ode := diff( y(x), x ) = y(x); 
icA := [0, 1]; 
icB := [0. 2]; 
icC := [0,-1]; 

DEplot( ode, y(x), x=0..2, [icA,icB,icC], arrows=slim, linecolor=blue, title="#23 (Section 9.1)" ); 
Mathematica : 

To plot vector fields, you must begin by loading a graphics package. 

«Graphics'PlotField' 

To control lengths and appearance of vectors, select the Help browser, type PlotVectorField and select Go. 
Clear[x, y, f] 
yprime = y (2 — y); 

pv = PlotVectorField[{ 1, yprime}, {x, —5, 5), {y, —4, 6], Axes — > True, AxesLabel — > {x, y}]; 

To draw solution curves with Mathematica, you must first solve the differential equation. This will be done with 
the DSolve command. The y[x] and x at the end of the command specify the dependent and independent variables. 
The command will not work unless the y in the differential equation is referenced as y[x]. 
equation = y'[x] == y[x] (2 — y[x]) ; 
initcond = y[a] == b; 

sols = DSolve[{equation, initcond), y[x], x] 

vals= {{0, 1/2}, {0,3/2}, {0,2}, {0,3}} 

f[{a_, b_}] = sols[[l, 1,2]]; 

solnset = Map[f, vals] 

ps = Plot [Evaluate [solnset, {x, —5, 5}]; 

Show[pv, ps, PlotRange — > {—4,6}]; 

The code for problems such as 31 & 32 is similar for the direction field, but the analytical solutions involve 
complicated inverse functions, so the numerical solver NDSolve is used. Note that a domain interval is 
specified. 

equation = y'[xj == Cos[2x — y[x]] ; 
initcond = y [0] == 2; 

sol = NDSolve [{equation, initcond}, y[x], {x, 0, 5}] 
ps = Plot[Evaluate[y[x]/.sol, {x, 0, 5}]; 

N[y[x] /. sol/.x -► 2] 

Show[pv, ps, PlotRange — > {0,5}]; 

Solutions for 34 can be found one at a time and plots named and shown together. No direction fields here. 

For 34, the direction field code is similar, but the solution is found implicitly using integrations. The plot 
requires loading another special graphics package. 

«Graphics' ImplicitPlot' 

Clear[x,y] 

solution[c_] = Integrated (y — 1), y] == Integrate[3x 2 + 4x + 2, x] + c 

values = {—6, —4, —2, 0, 2, 4, 6}; 

solns = Map[solution, values]; 

ps = ImplicitPlot[solns, {x, —3, 3}, {y, —3, 3}] 

Show[pv, ps] 
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27. 28. 




35. ^ = 2xe x \y(0) = 2 => y n+ i = y n + 2x n e x »dx = y n + 2x n e x »(0.1) = y n + 0.2x n e x « 

On a TI-84 calculator home screen, type the following commands: 

2 STO > y: 0 STO > x: y (enter) 

y + 0.2*x*e A (x A 2) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last value displayed gives yEuier(l) ~ 3.45835 

The exact solution: dy = 2xe x2 dx => y = e x2 + C; y(0) — 2 = e H + C => C — 1 => y = 1 + e x ' 

=> y exact (1) = 1 + e « 3.71828 

36. 37 = 2y 2 (x - 1), y(2) = ~\=> y n +i = y n + 2y 2 (x n - l)dx = y n + 0.2y 2 (x n - 1) 

On a TI-84 calculator home screen, type the following commands: 

—0.5 STO > y: 2 STO > x: y (enter) 

y + 0.2*y 2 (x — 1) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last value displayed gives yEuier(2) « —0.19285 

The exact solution: ^ = 2y 2 (x — 1) => p- = (2x — 2)dx => — f=x 2 — 2 x + C => 3 = _ x 2 + 2x + C 
y(2) = 4 =► 3 T /2 = -(2) 2 + 2(2) + C = C => C = —2 => 3 = —x 2 + 2x — 2 => y = _ x2+ ‘ 2x _ 2 

= — ( 3 ) : + 2 ( 3 ) — 2 = ”‘ 0 ' 2 
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37 - s = #’ y > °> y(°) = 1 =* y n +i = y n + #dx = Yn + ^(o.i) = Yn + o.i^ 

On a TI-84 calculator home screen, type the following commands: 

1 STO > Y : 0 STO > x: y (enter) 

y + 0.1*(v^ 7 y) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last value displayed gives Y E U ier(l) ~ 1.5000 

The exact solution: dy = ^dx => y dy = ^/x dx \ = |x 3/2 + C; = y = \ = |(0) 3/2 + C => C = \ 

=► i = |x 3/2 + y = yf fx 3 ' 2 + 1 =► y exact (l) = yj\{ l) 3/2 + l w 1.5275 

38. | = 1+ y 2 , y(0) = 0 => y n+1 = y n + (1 + y 2 )dx = y„ + (1 + y 2 )(0.1) = y n + 0.1(1 + y 2 ) 

On a Tl-84 calculator home screen, type the following commands: 

0 STO > y: 0 STO > x: y (enter) 

y + 0. 1 *( 1 + y 2 ) STO > y: x + 0.1 STO > x: y (enter, 10 times) 

The last value displayed gives Y Euier(l) ~ 1.3964 

The exact solution: dy = (1 + y 2 )dx => } , = dx => tan _1 y = x + C; tan _1 y(0) = tan^'0 = 0 = 0 + C=>C = 0 

=> tan~'y = x =>• y = tan x => y ex act(l) = tan 1 w 1.5574 

39. Example CAS commands: 

Maple : 

ode := diff( y(x), x ) = x + y(x);ic := y(0)=-7/10; 
xO := -4;xl := 4;y0 := -4; yl := 4; 
b := 1; 

PI := DEplot( ode, y(x), x=x0..xl, y=y0..yl, arrows=thin, title="#39(a) (Section 9.1)" ): 

PI; 

Ygen := unapply( rhs(dsolve( ode, y(x) )), x,_Cl ); # (b) 

P2 := seq( plot) Ygen(x.c), x=x0..xl, y=y0..yl, color=blue ), c=-2..2 ): # (c) 

display! [P1,P2], title="#39(c) (Section 9.1)" ); 

CC := solve! Ygen(0,C)=rhs(ic), C ); # (d) 

Ypart := Ygen(x.CC); 

P3 := plot( Ypart, x=0..b, title="#39(d) (Section 9.1)" ): 

P3; 

euler4 := dsolve( {ode,ic}, numeric, method=classical[foreuler], stepsize=(xl-xO)/4 ): # (e) 

P4 := odeplot( euler4, [x,y(x)], x=0..b, numpoints=4, color=blue ): 
display( [P3,P4], title="#39(e) (Section 9.1)" ); 

euler8 := dsolve( {ode,ic}, numeric, method=classical[foreuler], stepsize=(xl-x0)/8 ): # (f) 

P5 := odeplot( euler8, [x,y(x)], x=0..b, numpoints=8, color=green ): 

eulerl6 := dsolve( {ode.ic}, numeric, method=classical[foreuler], stepsize=(xl-x0)/16 ): 

P6 := odeplot( eulerl6, [x,y(x)], x=0..b, numpoints=16, color=pink ): 

euler32 := dsolve( {ode.ic}, numeric, method=classical[foreuler], stepsize=(xl-xO)/32 ): 

P7 := odeplot( euler32, [x,y(x)j, x=0..b, numpoints=32, color=cyan ): 
display( [P3,P4,P5,P6,P7], title="#39(f) (Section 9.1)" ); 

< < N | h | 'percent error' >, # (g) 

< 4 | (xl-x0)/ 4 | evalf[5]( abs(l-eval(y(x),euler4(b))/eval(Ypart,x=b))*100 ) >, 

< 8 | (xl-x0)/ 8 | evalf[5]( abs(l-eval(y(x),euler8(b))/eval(Ypart,x=b))*100 ) >, 

< 16 | (xl-x0)/16 | evalf[5]( abs(l-eval(y(x),eulerl6(b))/eval(Ypart,x=b))*100 ) >, 

< 32 | (xl-x0)/32 | evalf[5]( abs(l-eval(y(x),euler32(b))/eval(Ypart,x=b))*100 ) > >; 
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39-42. Example CAS commands: 

Mathematica : (assigned functions, step sizes, and values for initial conditions may vary) 

Problems 39 - 42 involve use of code from Problems 23 - 34 together with the above code for Euler's method. 



9.2 FIRST-ORDER LINEAR DIFFERENTIAL EQUATIONS 

1. x|+ y = e ^| + (I)y=f,P(x ) =I,Q( x) =f 

f P(x) dx = Ji dx = In | x | = In x, x > 0 => v(x) = e/ p(x)dx = e lnx = x 
y=v|)/^ x )QWdx=i/x(f)dx=i(e x + C) = ^,x >0 

2 . e x g + 2 e x y = 1 => g + 2 y = e" x , P(x) = 2 , Q(x) = e~ x 

J P(x) dx = J 2 dx = 2 x v(x) = e i P(x) dx = e 2x 

y = /e 2x • e x dx = ^ J e x dx = ^ (e x + C) = e~ x + Ce- 2x 

3 . xy ' + 3 y = 5 |x , x > 0 ^ g + (f ) y = ^ , P(x) = 3 , Q(x) = 

J 5 dx = 3 In |x| = In x 3 , x > 0 => v(x) = e lnxi = x 3 

y = fx 3 (^) dx = i / sin x dx = i (_ cos x + C) = , x > 0 

4 . y' + (tan x) y = cos 2 x, — |<x<|=^$[-|- (tan x) y = cos 2 x, P(x) = tan x, Q(x) = cos 2 x 
J\an x dx = J dx = — In |cos x| = In (cos x) _1 , — f < x < f => v ( x ) = e ln(cosx) = (cos x )^ 1 
y = ^0^-1 f (cos x )- 1 • cos 2 x dx = (cos x) J cos x dx = (cos x)(sin x + C) = sin x cos x + C cos x 



dy 

dx 



2y=l — ^,x>0=> 



+ (i) y = x - ?■* P(x > = x > QW = x 



J= dx = 2 In | x | = In x 2 , x > 0 => v(x) = e lnx2 = x 2 



y = 



r / x2 (I-i)d X =i/( x -l)dx=i(f- x + C ) = !-! + §, x>0 



6 . 



(l+x)y' + y=^^| + ( T ^)y 



jA 

1+ x 



,P(x) 



1 

1+x 



, Q(X) = 



_A_ 

l+x 



/ 



y = 



dx = In ( 1 + x), since x > 0 



dx = 

1 +x 



l+x 



/ <> + *> ( 4 ) 



>■ v(x) = e ln(1+x) = 1 

fV* d X =( r ^)(|x 3 / 2 + C) 



2x 3 / 2 , C 

3(1 +x) f l+x 



7. S - \ y = l e x/2 =» P(x) = - Q(X) = \ e x/2 =>• f P( X ) dx 

^ y = ^ J e * /2 (i e x/ 2 ) dx = e x / 2 f I dx = e x/2 (| x + C) 



-|x v(x) = e x/2 
| xe x/2 + Ce x/2 



8. g + 2y = 2xe- 2x => P(x) = 2, Q(x) = 2xe- 2x => / P(x) dx = Jl dx = 2x =>• v(x) = e 2x 

=$■ y = f s 2x (2xe~ 2x ) dx = ^ fix dx = e~ 2x (x 2 + C) = x 2 e~ 2x + Ce 2x 

9. g - (i) y = 2 In x =*► P(x) = - \ , Q(x) = 2 In x => f P(x) dx = - f ± dx = - In x, x > 0 

=> v(x) = e~ lnx = i => y = xj'(l) (2 In x) dx = x [(In x) 2 + C] = x(ln x) 2 + Cx 
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t°. g + (?) y = «S|* , x > 0 => P(x) = l , Q(x) = «§? => f P(x) dx = fl dx = 2 In |x| = In x 2 , x > 0 

=> v(x) = e lnx2 = x 2 => y = 4 Jx 2 (^) dx = 4 /cos x dx = ^ (sin x + C) = + c 

H. $ + (r^r) s = (T^TF =* P(t) = rrr , Q(t) - (r^TF =* / P(t) dt = f± dt = 4 In |t - 1| = ln(t- l) 4 

=* v(t) = e ln(t ~ 1)4 = (t - l) 4 =► s = ^ f(t-l) 4 [^] dt = ^ /( t 2 - 1) dt 

1 ( t c 

= (t-1) 4 ^ 3 ~ 1 = 3(t— I) 4 _ (t- l) 4 ' ft - l) 4 

12. (1+1)1+28 = 3(1+1)+^ =s- $ + (t+r) 8 = 3 + (hT? =}► P(t)= ^,Q(t) = 3 + (t+lr 3 

=> / P(t) dt = f dt = 2 In |t + 1| = ln(t + l) 2 => v(t) = e ln(t+1 > 2 = (t + l) 2 

=> s = /(t + l) 2 [3 + (t + 1) 3 ] dt = J [3(t + l) 2 + (t + 1) 3 ] dt 

= (T+TF [9 + l) 3 + In |t + 1 1 + C] = (t + 1) + (t + 1) 2 ln(t + 1) + )2 , t > — 1 

13. ^ + (cot 0) r = sec 0 =F P(0) = cot 0, Q(0) = sec 0 =F / P(0) d0 = /cot 0 dO = In |sin 0j => v(0) = e ln|si “ fl| 

= sin 0 because 0 < 0 < | =>■ r = — L / ( s i n g)( sec 0 ) d 0 = — L / j an g ^g — _L_ q n | sec $| c) 

= (esc 9 ) (In |sec 9\ + C) 

14. tan 9 (§ + r = sin 2 0 => ^ + ^ = ^| =>• ^ + (cot 0)r = sin 9 cos 9 => P(0) = cot 0, Q(0) = sin 0 cos 0 

=>• J P(0) d 9 = f cot 0 d0 = In |sin 0| = In (sin 0) since 0 < 0 < | =>■ v(0) = e ln(sm63 = sin 0 

=» r = site f( S4n V 0 cos e (iib) + C) = + £s 

15. | + 2y = 3 =>• P(t) = 2, Q(t) = 3 => f P(t) dt = / 2 dt = 2t => v(t) = e 2t =>• y = ^ / 3e 2t dt 

~ e® (I e2t + C) ; y(0) = 1 ^!+C=l=»C = -i=> y = § - § e^ 

16. | ^ = t 2 P(t) = \ , Q(t) = t 2 => / P(t) dt = 2 In |t| => v(t) = e ln ' 2 = t 2 =» y = ^ /( t 2 ) (t 2 ) dt 

= F f t4 dt = F (l + C ) = T + F ; y( 2 ) =1 => f + 4 = 1 C = _ T y=F _ 5F 

17 . | + Q) Y = ^ => W) = | , Q(0) = ^ =+ / P(0 d0 = In |0| => v(0) = e 1 "' 8 ' = |0| 

^y^ /|0| d0= i/0(^) d0for0^O =* y = i/sin 0 d0 = i (- cos 0 + C) 

= - i cos 0 + f ; y (|) = 1 =+ C = | ^y = -l C os0+^ 

18. — (|) y = 0 2 sec 0 tan 0 => P(0) = — | , Q(0) = 0 2 sec 0 tan 0 => J' P(0) d0 = —2 In |0| => v(0) = e- 21 ” 1 " 1 
= 0 -2 =>• y = glj J" (0 -2 ) (0 2 sec 0 tan 0) d0 = 0 2 J* sec 0 tan 0 d0 = 0 2 (sec 0 + C) = 0 2 sec 0 + C0 2 ; 

y (f ) = 2 => 2 = (f ) (2) + C (f ) => C = if - 2 =*► y = 0 2 sec 0 + (« - 2) 0 2 



=» jy _ 2 f x(x+ !> ] __ 

+ x ^y- x+i ^ dx z |_x+i jy- 

f P(x) dx = f — 2x dx = -x 2 => v(: 
= e x2 [^+c]=-^+Ce x2 ;; 



e x2 dy /-* 

y = (F+tf =► dF - 2x y : 
v(x) = e- x2 => y= ^ 
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20. g + xy = x =>• P(x) = x, Q(x) = x => f P(x) dx = fx dx = f => v(x) = e x2 / 2 =>• y = \ f e x2/2 • x dx 

= ^ (e xV2 + C) = 1 + ;y(0) = -6 ^ 1 + C = -6 ^ C = -7 ^ y = 1 - ^ 

21. f t - ky = 0 =x P(t) = -k, Q(t) = 0 => J P(t) dt = § k dt = -kt => v(t) = e~ kt 

=> y = f (e” kt ) (0) dt = e kt (0 + C) = Ce kt ; y(0) = y 0 =xC = y 0 =xy = y 0 e kt 

22. (a) | + £ u = 0 =» P(t) = | , Q(t) = 0 => f P(t) dt = / £ dt = £ t = £ =s- u(t) = e k '/“ 

=> y = ^f$ /e kt/m • 0 dt = jL, ; u(0) = u 0 => cE ^ 7il = u 0 => C = u 0 => u = u 0 

(b) ^ = — - u=>- — = — -dt=>ln U = — - t + C=» U = e _(k / m ) t+c =>■ u = • e c . Let e c = Ci . 

v/ dtmum m 1 

Then u = ^ • Ci and u(0) = u 0 = • Ci = Ci. So u = u 0 

23. x J' ^ dx = x (In |x| + C) = x In |x| + Cx => (b) is correct 



24. 



1 

COS X 



I cos x dx = — — (sin x + C) = tan x H — — 

J COS X v J COS X 



=> (b) is correct 



25. Steady State = ^ and we want i = | (g) =>• | (^) = ^ (1 — e Rt/L ) =X \ — 1 — e Rt/L =X — \ = — e Rt / L 

=> In | ^ => — g In \ = t => t=gln2sec 

26. (a) f + £ i = 0 =X \ di = - | dt =X In i = - f + Cj =X i = e c 'e- Rt T = Ce- Rt / L ; i(0) = I =>■ 1 = C 

=> i = Ie~ Rt/L amp 

(b) |l = I e Rl/L =>■ e~ Rt,/L = 1 => — ^ = In i = — In 2 => t = h l n 2 sec 

(c) t = h => i = I e (-Rt/L)(L/R) _ j e -t am p 

27. (a) t — 3 j k =$■ i = g (1 — e( _R/L) ( 3L / R )) = L (1 — e -3) ~ 0.9502 ^ amp, or about 95% of the steady state value 

(b) t = ^ => i = X (i _ e (— r/l)( 2 l/r)^ _ v ^ _ e - 2 ^ py 0.8647 ^ amp, or about 86% of the steady state value 



28. (a) 



(c) 



di 

dt 



!i = 



P(t) = l , Q(t) = l 



^ dt = — 

L L L 



1 = 



1 

e Rt/L 



f ?«)“' = / 

J> L (r) * = * [I = ,,,L (?) + c] = l + Ce-«™ 



v(t) = e 1 



,Rt/L 



(b) i(0) = 0 ^ £+C = 0^C = -£ => i = 



i _ V _ V p— Rt/L 



R 



R 



i=S 



gj = 0 => g| + R i = 0 + ( R ) (g) = ^ => i = ^ is a solution of Eq. (1 1); i = Ce (R/L )‘ 



29. y' — y = — y 2 ; we have n = 2, so let u = y 1 2 — y *. Then y = u 1 and ^ = — ly 2 gj => gj = — y 2 ^ 
=> -u" 2 ^ - u- 1 = -u" 2 ^ + u = 1. With e/ dx = e x as the integrating factor, we have 

eX (s + u ) = ®( eXu ) = eX - Integrating, we get e x u = e x + C =>• u = 1 + § = j => y = 



30. 



y' - y = xy 2 ; we have n = 2, so let u = y” 1 . Then y = u^ 1 and | = -y^ => g = -y 2 ^ = -u" 2 ^. 
Substituting:— u _2 jy — u~' = xu -2 => ^ + u = — x. Using ef dx = e x as an integrating factor: 



[ (g +u) = ^(e x u) = — x e x =^e x u = e x (l - x) + C 



=)> u = 






- xe x + C 



31. xy ' + y = y 

du _ o 2 dy 
dx y dx 



- 2 => y' + (i)y = (I)y -2 . Let U = yH^ = y 3 => y = u 1/3 

=*-y' = S = d)(S(r 2 ) = d)(|)(u-”). THUS we have 



and 



2 



= U 



- 2/3 
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(5) (gj) ( u 2/3 ) + (D u ‘ /3 = (x) u 2/3 => S + (x) u = (x) h The integrating factor, v(x), is 
eix dx = e 31nx = e lnx3 = x 3 . Thus ( g:(x 3 u) = (f)x 3 = 3x 2 => x 3 u = x 3 +C=>u=l + ^=y 3 

=>y=(i + §) 1/3 

32. x 2 y' + 2xy = y 3 ^ y' + (|)y = (L)y 3 . P(x) = (jj),Q(x) = (4), n = 3. Let u = y 1 - 3 = y" 2 . 

Substituting gives ^ + (— 2)( = )u = — 2(^) => ^ + (^)u = -jr. Let the integrating factor, v(x), be 
e/(4) dx = e lnx 4 = x -4 . Thus (^(x -4 u) = — 2x -6 => x" 4 u = jX -5 + C => u = L|+Cx 4 = y -2 
^y=(li+cx 4 )- 1/2 

9.3 APPLICATIONS 

1. Note that the total mass is 66 + 7 = 73 kg, therefore, v = voe _ ^ m ^ => v = 

(a) s(t) = / 9e _3 - 9t / 73 dt = _2i|0 e -3.9t/73 + c 

Since s(0) = 0 we have C = 422 anc j li m s(t) = lim 4P (l — e _3 - 9t / 73 ) 

Woo Woo v ' 

The cyclist will coast about 168.5 meters. 

(b) 1 = 9e~ 3 - 9t / 73 =s> L|l = in 9 =>. t = ^LL 9 « 41.13 sec 
It will take about 41.13 seconds. 

2 V = v 0 e~( k,/m)t => V = 9 e — (59-000/51, 000,000)t v _ g e -59t/51,000 

(a) s(t) = / 9e — 59t / 5 1 000 dt = e -59t/5i.ooo + c 

Since s(0) = 0 we have C = and lims(t) = lim ^“992 (1 - e -59t/5i,ooo\ = iswwoo _ 7780 m 

Woo Woo v ' 

The ship will coast about 7780 m, or 7.78 km. 

(b) 1 = 9e~ 59t/51 ’°°° => = i n 9 ^ t = 51 ’°°° ln9 w 1899.3 sec 

It will take about 31.65 minutes. 

3. The total distance traveled = ^ =>■ ^ 2 ' 75) [ 39,92) = 4.91 => k = 22.36. Therefore, the distance traveled is given by the 

function s(t) = 4.91 (l — e^ 22 - 36 ^ 9 - 92 ) 1 ) . The graph shows s(t) and the data points. 




9e -3.9t/73 

= » 168.5 



4. ^ = coasting distance => (°- 8 ( U 49 - 9Q ) = 1.32 =>■ k = ^ 

We know that = 1.32 and L = = §2. 

Using Equation 3, we have: s(t) = ^ (l — e~! k/n 9 l ) — 1 .32 ( 1 — e _20t/33 ) « 1.32(1 — e _a606t ) 
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5. y = mx => | = m => xy x2 y = 0 => y' = So for 
orthogonals: ^ ^ => y dy = — x dx => + ^ = C 

=> x 2 + y 2 = Ci 



6. y = cx 2 =>■ = c =>■ x y x 4 2xy = 0 => x 2 y' = 2xy 

=> y' = y • So for the orthogonals: gj = — ^ 

=> 2ydy = -xdx =^y 2 = -^+C=>y = ± ^ f + C, 

C > 0 



7. kx 2 + y 2 = 1 => 1 - y 2 = kx 2 => = k 

=> x2(2y)y ';i 1 - y2)2x = 0 =► — 2yx 2 y' = (1 - y 2 )(2x) 
=> y' = ^ _o x |, 2 "' x) = ^_ X y - . So for the orthogonals: 

S = i ^^ 11 7 fld J = xdx ^ 1 "y-? = T+ c 



8. 2x 2 + y 2 = c 2 =» 4x + 2yy' = 0 => y' = - ^ — . For 

orthogonals: ^ ^ — = ^ => In y = |ln x + C 

=> In y = In x 1/2 + InCi => y = Ci |x| 1/2 



n —x v e V -y(e x )(— l) n 

9. y = ce x => = c => — y /V, = 0 

J e (e- x ) 



=> e x y' — —ye x => y ' — — y. So for the orthogonals: 
| = i =>ydy = dx^ ^ =x + C 
=> y 2 = 2x + Ci => y = d= yj 2x + Ci 





y 




y 
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10. y = e kx =$■ In y = kx =>■ = k => ^ x 2 = 0 

=> ^ j y' — In y = 0 => y' = So f or the orthogonals 

S = ^=^y ln y d y = - xdx 

=^b 2ln y-?(y 2 ) = (-5 x2 )+ c 

=S> y 2 ln y - ^ = -x 2 + Ci 

11. 2x 2 + 3y 2 = 5 and y 2 = x 3 intersect at (1, 1). Also, 2x 2 + 3y 2 = 5 => 4x + 6y y' = 0 => y' = — =>■ y'(l, 1) = — | 

y 2 = x 3 => 2yi y [ = 3x 2 =>■ y{ = => y [ ( 1 , 1) = |. Since y' ■ y[ = (— |) (|) = — 1, the curves are orthogonal. 

2 2 

12. (a) x dx + y dy = 0 =>■ y + \ = C is the general equation 

of the family with slope y' = — For the orthogonals: 
y / = y ^ dy = ^ ^ In y = in x + C or y = Cjx 
(where Ci = e c ) is the general equation of the 
orthogonals. 



(b) x dy - 2y dx = 0 =>• 2y dx = x dy => ^ = f 

=> i = f =>• ±ln y = In x + C =>■ y = Qx 2 is 

the equation for the solution family. 

it In y — lnx = C=>i — — - =0^>y'= — 

2 J 2 y x J x 

=> slope of orthogonals is ^ = — j- 
=> 2y dy = — x dx => y 2 = — y + C is the general 
equation of the orthogonals. 



13. Let y(t) = the amount of salt in the container and V(t) = the total volume of liquid in the tank at time t. 
Then, the departure rate is ^4 (the outflow rate). 

(a) Rate entering = = 10 lb/min 

(b) Volume = V(t) = 100 gal + (5t gal - 4t gal) = (100 + t) gal 

(c) The volume at time t is (100 + 1) gal. The amount of salt in the tank at time t is y lbs. So the 
concentration at any time t is 10 ^ lbs/gal. Then, the rate leaving = 10 ^ (lbs/gal) • 4 (gal/min) 

= m+i lbs/min 

W ! = 10 -Iofe =* f + (iW+t)y= 1 0 =► P(t)=To^,Q(t)=10 =* /P(t) dt = /to^r dt 
= 4 In (100 + t) => v(t) = e 41n(100+t > = (100 + t) 4 =► y = <10^1 /(100 + 1) 4 (10 dt) 

= (TMT1F + c) = 2( 1 00 + 1 ) + ; y(0) = 50 =* 2(100 + 0) + = 50 

=► C = -(150X100) 4 =► y = 2(100 + 1) - %°> ( l°y =* y = 2(100 + 1) - -^0 
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(e) y(25) = 2(100 + 25) - « 188.56 lbs =* concentration = v y j^ c » ^ » 1.5 lb/gal 

14. (a) ^ = (5 - 3) = 2 => V = 100 + 2t 

The tank is full when V = 200 = 100 + 2t => t = 50 min 

(b) Let y(t) be the amount of concentrate in the tank at time t. 

dy _ ( 1 Jb (r gaA _ ( y lb \ Ggal\=>dy = 5_ 3 / y \ dy , 3 _ 5 

dt y 2 gal J y min J yl00 + 2tgal J y min J dt 2 2 \50 + t/ dt 2(t + 50)' 2 

Q(t) = P(t) = I (joVi) =>■ / P(t) dt = dt = fin (t + 50) since t + 50 > 0 

v(t) = e i' P(t,dt = eJ ln(t+5 °) = (t + 50) 3/2 

y(t) = (t ^ 50 p /i(t + 5 °) 3/2 dt = (t + 50) - 3/2 [ (t + 50) 5/2 + c] =► y(t) = t + 50 + 

Apply the initial condition (i.e., distilled water in the tank at t = 0): 

y(0) = 0 = 50 + => C = — 50 s/2 => y(t) = t + 50 — - ^ . When the tank is full at t = 50, 

y(50) = 100 — pA ~ 83.22 pounds of concentrate. 



15. Let y be the amount of fertilizer in the tank at time t. Then rate entering — I - I = 1 



and the 



volume in the tank at time t is V(t) =100 (gal) + [1 (gal/min) — 3 (gal/min)]t min = (100 — 2t) gal. Hence 



rate out = I 

=► P(t) = 



y ) 

100-2t ) 



3 = TO0=2t lbs/mil1 



dy 



dt ( ^ 100 — 2 tl 



3y 



lbs/min 



too— 2 t 
= (100 - 2t) -3 / 2 



,Q(t) = l =► f P(t)dt= / 



3 Hf . _ 3 ln(100— 2t) 
100 — 2 t UL -2 



* id + (too^s) y - 1 

v(t) = e (-31n(100-2t))2 



f (100 - 2t) -3 / 2 dt = (100 - 2t) 



3/2 



Let dy _ q dy _ _ 2 
dt — u ^ dt — z 



y (too -2tr 3 / 2 

= (100 - 2t) + C(100 - 2t) 3 / 2 ; y(0) = 0 =+ L100 - 2(0)] + C[100 - 2(0)] 3 / 2 
=> C= -(100 r 1/2 = - i => y = (100 - 2t) - (10Q 7 0 2t)3/2 
= -2 + 3v/l °g~ 2t = 0 =+ 20 = 3 1/ 100 — 2t => 400 = 9(100 - 2t) 

=> t « 27.8 min, the time to reach the maximum. The maximum amount is then 
y(27.8) = [100 - 2(27.8)] - t 100 - 2 ^ 7 -^ 2 ~ 14.8 lb 



- 2(100 — 2 t )~ 1 2 



C(100) 3 / 2 = -100 



400 = 900 - 18t 



(]) (100 — 2 t) 1 / 2 (— 2 ) 

10 

-500 = -18t 



16. Let y = y(t) be the amount of carbon monoxide (CO) in the room at time t. The amount of CO entering the 
room is (^55 x TT)) = TM) ft 3 /min, and the amount of CO leaving the room is (pp) (y/j) = tTooo ft 3 /inin. 



Thus, $ = 12 



dy 



1000 15,000 



+ 



12 



P(t) = 



, Q(t) = 



12 

1000 



15,000 

v(t) = e 1715 ’ 000 



15,000 y — 1000 * w — 15,000 

=*► y = fj^o e 1 / 15 - 000 dt => y = eW 15 >°°° (l 2 ™ e t/i5,ooo + c ) = e -t/i 5 ,ooo (180e ./u,ooo + c) . 

y(0) = 0 =+ 0=1(180 + 0 => C = — 180 => y = 180 - 180e- ,yi5 ’ 000 . When the concentration of CO is 0.01% 
in the room, the amount of CO satisfies 4 y J(J = => y = 0.45 ft 3 . When the room contains this amount we 

have 0.45 = 180 - lSOe^ 15 ’ 000 => ± 2+55 = e -t/i5,ooo ^ t = -15,000 In « 37.55 min. 



180 



9.4 GRAPHICAL SOUTIONS OF AUTONOMOUS EQUATIONS 



L y' = (y + 2)(y-3) 

(a) y = — 2 is a stable equilibrium value and y = 3 is an unstable equilibrium. 

(b) y" = (2y - l)y' = 2(y + 2) (y - \) (y - 3) 



y’> 0 


i y’< 0 


i y’> 0 


i ^ 


a 4 1 i 


X 1 -w . 


-4 


+ * 0 1 2 


T 4 > 5 


y"< 0 


j y*>0 [ y'<0 


! y">0 



0.5 
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(c) 




2. y' = (y + 2)(y — 2) 

(a) y = — 2 is a stable equilibrium value and y = 2 is an unstable equilibrium. 

(b) y" = 2yy' = 2(y + 2)y(y — 2) 



y' > 0 


y' < 0 


y' > 0 


-3 * - 


[2 -1 1 


i° 1 \ 


3 * 


•< 

A 

o 


V 

o 


A 

o 


y" > 0 



(c) 




3. y' = y 3 -y = (y+l)y(y-l) 

(a) y = — 1 and y = 1 is an unstable equilibrium and y = 0 is a stable equilibrium value. 

(b) y" = (3y 2 - l)y' = 3(y + 1) (y + ^)y(y - (y - 1) 



y' < 0 


y' > 0 


y' < 0 


y' > 0 


-1 .5 * J 

y" < 0 


1 -0.5 ' 

y" > 0 


y" < ( 


y" > C 


O 

V 

io ^ 
r O . 


1.5 

y" > 0 



1 1 

a /3 y/Z 



y 
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y' = y(y-2) 

(a) y = 0 is a stable equilibrium value and y = 2 is an unstable equilibrium. 

(b) y" = (2y — 2)y' = 2y(y — l)(y — 2) 



y' > 0 


y' < 0 

. 


o 

A 


-l * | 


5 il 


i i 


-3 


A 

o 


V 

o 


O 

V 

v» 


V 

o 



y 



y> 0, y*> 0 
y'< 0. y* < 0 



y<o, y > o 



/ > 0, y <0 




y'= y^-y >0 



(a) There are no equilibrium values. 




y' = y - yy, y > 0 

(a) y = 1 is an unstable equilibrium. 



< b > y " = (> - A) y ' = (‘ - A)< y - AO = (A - IMA - >) 



y' < 0 


<< 

A 

o 


y' >0 






0.5 i 


1.5 


’ 2 


O 

A 


A 

o 


o 

A 

V> 




0.25 
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(c) 



y 




7. y' = (y- l)(y-2)(y-3) 

(a) y = 1 and y — 3 is an unstable equilibrium and y = 2 is a stable equilibrium value. 

(b) y" = (3y 2 - 12y + 1 1 ) (y - l)(y - 2)(y - 3) = 3(y - 1) (y - (y - 2) (y - (y 



1 


y<o 


• y>o i 


y'< 0 


y>o 


i 


0 




1? ! 2? 


; 3? 




4 




y”< 0 


1 y">0 1 y*<0 1 


y"> 0 i y"<0 1 


y"> 0 








6 -^=1.42 


6 + ^~2.58 







3 3 



(c) 



y 



4 

3.5 






— ■ — y'< 0, y"< 0 




y'<0,y">0 




— ■ — U y’> 0, > ,r < 0 




— y'>0,y'>0 


0.5 


N,^ y'<0,y ff <0 


-1 


1 2 3 ^ 


3 2 

= y - y 


(N 

II 



(a) y = 0 and y = 1 is an unstable equilibrium. 

(b) y" = (3y 2 - 2y)(y 3 - y 2 ) = y 3 (3y - 2)(y - 1) 



• VJ_ 

A 

o 


y' < 0 


y' > 0 


o 

V 

->> 


O 

A 

V» 


y" < 0 


t 


O 

A 

V> 



2/3 



(c) 



y 
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